Ground state entanglement of the BCS model 
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The concept of local concurrence is used to quantify the entanglement between a single qubit and 
the remainder of a multi-qubit system. For the ground state of the BCS model in the thermodynamic 
limit the set of local concurrences completely describe the entanglement. As a measure for the 
entanglement of the full system we investigate the Average Local Concurrence (ALC). We find that 
the ALC satisfies a simple relation with the order parameter. We then show that for finite systems 
with fixed particle number, a relation between the ALC and the condensation energy exposes a 
threshold coupling. Below the threshold, entanglement measures besides the ALC are significant. 
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Quantification of entanglement remains as a major 
challenge in quantum information theory. For a state of 
a bi-partite system, it is well known that a measure for 
the entanglement between the two subsystems is given by 
the von Neumann entropy y|. For a two-qubit system, 
the concurrence has been introduced as an alternative 
measure which is related to the von Neumann entropy in 
a bijective manner. In the language of spin-1/2 particles, 
concurrence can be described in terms of a time-reversal 
operation. Using this concept, generalisations of concur- 
rence have been proposed for multi-qubit systems and 
in particular they have been applied to quantify the en- 
tanglement of the ground state of the BCS model in the 
thermodynamic limit 0, Q , which breaks time-reversal 
symmetry due to broken gauge symmetry. 

In this Letter we adopt a different approach to inves- 
tigate the ground state entanglement of the BCS model. 
We will use the notion of local concurrence which is de- 
fined in analogy with the functional relation that exists 
between concurrence and von Neumann entropy for a two 
qubit system (cf. [f|). It is a measure of the entangle- 
ment between a single qubit and the remainder of the 
system. We then define an entanglement measure which 
is the Average Local Concurrence (ALC). The ALC satis- 
fies the properties of an entanglement monotone (EM): it 
vanishes for a state if and only if that state is a product 
state, is invariant under local unitary transformations, 
and does not increase on average under local operations 
assisted by classical communication. 

For multi-qubit systems the number of independent 
EMs is the same as the number of non-local invariants 
, growing exponentially with system size (2 L+1 — 3L — 2 
where L is the number of qubits 0, @)- It is therefore 
useful to identify EMs which can be related to physical 
aspects of the system. In the thermodynamic limit we 
will show that the ALC for the ground state of the BCS 
model displays a simple relationship with the magnitude 
of the order parameter. For finite systems the order pa- 
rameter vanishes, so we may use the ALC in lieu of the 
order parameter. Our investigation of a relationship be- 
tween the ground state ALC and the condensation energy 
exposes a threshold coupling which signifies the onset of 
entanglement not measured by the local concurrences. 



To study the entanglement of any pure state we can 
examine the von Neumann entropy. Consider a general 
quantum system comprised of L qubits and a partition 
into two subsystems denoted A and B. For any pure 
state density matrix p the entanglement Bab ifi) between 
A and B is given by the von Neumann entropy 

£ab(p) = -tr(p A log p A ) = -tr(p B logps) 

where the logarithm is taken base 2 and pA is the re- 
duced density matrix obtained from p by taking the par- 
tial trace over the state space of subsystem B. The re- 
duced density matrix pb is defined analogously. 

Now we make precise the definition of local concur- 
rence in a general context. Hereafter we will only deal 
with the particular case when the subsystem A denotes 
a single qubit (say the jth qubit) and B denotes the re- 
mainder of the system. For this case we will write Pj for 
Pa and £j(p) for £ab(p)- In such an instance we have 



A, log A: 



where A^ denote the 
1 



£ s (p) = -Atlog(At) 

two eigenvalues of pj. Using the fact that tr(/9j) 
and that the eigenvalues of pj lie in the interval [0, 1] 
means that we can always parameterise them as A^ — 

(l ± Jl - C? ) /2 with Cj e [0, 1]. For a two-qubit sys- 
tem C — C\ — C2 is precisely the concurrence 2], so it 
is natural to call the Cj for the L-qubit system the lo- 
cal concurrences (cf. jjj). In terms of Pauli matrices it 

can be determined that pj = (i + J2 a =x y z { a f) a jJ 
giving the local concurrence as 



(i) 



For the ground state of the BCS model we will establish 
a correspondence between each local concurrence and a 
certain correlation function Cj (see © below) describing 
the fluctuation in the Cooper pair occupation numbers. 
An advantage of this approach is that it applies equally 
to the thermodynamic limit where gauge symmetry is 
broken, and for finite systems where there is no broken 
symmetry. We obtain analytic results for the ALC in 
two extreme cases; (i) the thermodynamic limit and (ii) 
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the case of a single Cooper pair in a system with two 
single particle energy levels. Between these extremes we 
investigate the ALC in terms of the exact solution of the 
model provided by the Bethe ansatz , which facilitates 
the calculation of correlation functions [Tol Hi) and in 
turn the ALC. 

The reduced BCS model. The reduced BCS Hamilto- 
nian has received much attention as a result of the effort 
to understand pairing correlations in nanoscale metallic 
systems [lj, UM ■ The Hamiltonian reads 



BCS 



L 



3 3 



-dx± 



C 3 + C 3- 



Ck-Ck+- 



(2) 



Above, j — 1,...,L labels a shell of doubly degenerate 
single particle energy levels with energies 6j , d is the mean 
level spacing and A is the dimensionless coupling. The 
operators Cj±, cj_|_ are the annihilation and creation op- 
erators for electrons at level j where the labels ± refer to 
pairs of time-reversed states and Uj — c^ + cj + + aj_Cj_ 
is the electron number operator for level j. Throughout 
we will only consider systems at half-filling. The physical 
properties predicted by the model are quite different in 
the superconducting (SC) regime (d <C A, where A is 
the bulk gap, defined below) and the fluctuation domi- 
nated (FD) regime (d ^> A), the latter being the case 
for nanoscale systems. For the SC regime the variational 
BCS ansatz using mean-field theory can be used to de- 
termine the ground state properties. However, the mean- 
field approximation is not justified in the FD regime, 
where quantum fluctuations are significant. 

The condensation energy E c is defined as the energy 
loss relative to the energy of the uncorrelated state (i.e. 
the ground state energy at A = 0, or equivalently the en- 
ergy expectation value of the Fermi sea). In the language 
of [13, it is equivalent to the entanglement gap for this 
model. Intuitively, the entanglement gap gives an indi- 
cation of the ground state entanglement of the system. 
By definition, it is zero if and only if the ground state 
is not entangled. It is thus desirable to determine how 
the entanglement gap (or equivalently the condensation 
energy) relates to EMs. It is known that the condensa- 
tion energy is extensive in the SC regime, intensive in 
the FD regime, but with entirely smooth behaviour in 
the crossover [2, 0] . Below we will use a relation be- 
tween condensation energy and the ALC to establish that 
a threshold coupling exists which marks qualitative dif- 
ferences in the ground state in terms of entanglement. 
An important quantity in our subsequent analysis will 
be played by the dimensionless condensation energy per 
electron, defined by E = E c /lodL. 

Next we discuss the decomposition of the Hilbert 
space into subsystems. At each energy level €j there 
are four independent states; |0) , c+ |0) , c_ |0) , c+cl |0). 
These states serve as a ququadrit, which can be further 
decomposed into two qubits through the identification 
10) = 100) , C++ |0) = |10) , cl |0) = |01) , c\J_ |0) = |11) 
[la Hq . In the ground state all electrons are paired, so 



there is zero probability of observing a single electron at 
any level tj . Thus for the ground state each level tj gives 

rise to a two-state system with basis |0) , c^cL |0), and 
each ququadrit serves as an effective qubit. 

The gran d canonical ensemble. The conventional BCS 
theory |l7| employs a grand canonical ensemble, where 
the electron number is not fixed, using the well-known 
variational ground state ansatz 



|BCS) 



L 



u 3 I 



(3) 



with 



real and satisfying u 2 + v 2 = 1. Including 



only those levels within the cut-off given by the Debye 
frequency ojd (i.e. < ljd where the Fermi level is 
ep = 0), minimisation of the expectation value of the 
energy for © gives 



A 2 /(e 2 + A 2 ) 



(4) 



where A = ui d / sinh(l/A) is the bulk gap. Unless 
stated otherwise we assume that the levels Cj are uni- 
formly distributed. It can then be deduced that E = 
(coth(l/A) - l)/2. 

With respect to the decomposition of J3J into 
ququadrits discussed above, it is a product state and thus 
not entangled. There is however entanglement within 
each ququadrit subsystem. Expressing the state of a 
ququadrit as u |00) +v e 10 |11) it is easily determined that 
the concurrence is C — 2\uv\. Thus for the state (J3J) we 
may define L local concurrences Cj associated with each 
level j, which quantify the entire entanglement content 
of the state. Remarkably, most of the 2 i+1 — 3L— 2 in- 
dependent EMs are zero. Note that the definition of lo- 
cal concurrence here is not the same as the definition of 
partial concurrence in based on the notion of broken 
time-reversal symmetry of © (see also 3]). A simple 
choice for an EM which reflects the overall ground state 
entanglement is the ALC, C = L^ 1 J2f=i Or We remark 
that C only quantifies bi-partite entanglement, and not 
multi-partite entanglement 0, Q, 13 • However (0 has no 
multi-partite entanglement. In the thermodynamic limit 
d — * with Ld = 2lod finite we can compute the ALC as 



- J" 
C = — 



f(eMe)de 



(5) 



with /(e) = A/(Ve 2 + A 2 ) and fi(e) > a density func- 
tion for the distribution of the single particle energy 
levels. For the case //(e) = 1 it is straightforward to 
determine that the ALC is C = l/(Asinh(l/A)). Re- 
calling j 1 ll ] that the order parameter is given by A = 



\dz 



we have |A| = XlujjC. This last rela- 



tion reflects the fact that superconducting order arises 
from the instability of the Fermi sea due to Cooper pair- 
ing, which also results in the emergence of entanglement 
in the ground state. Additionally we find 2dE/d\ = C . 



0.5 1 1.5 2 ^ 2.5 3 3.5 4 



0.5 1 1.5 2 ^ 2.5 3 3.5 4 



FIG. 1: The ALC as a function of the dimensionless coupling 
A in the thermodynamic limit for various distributions /x(e) 
of the single particle levels. The results shown are for, in 
increasing order of local density about the Fermi level, fi(e) = 
e 2 (dash), /x(e) = |e] (solid), fj,{e) = 1 (bold), /x(e) = uP D - e 2 
(dash-dot) and /x(e) = ojd — |e| (dot). For each of these cases 
the expression @ can be evaluated analytically. 

For the canonical case (i.e. a finite system with fixed elec- 
tron number) the order parameter vanishes, but one can 
alternatively study the extent to which 'superconducting 
order' survives through the study of certain correlation 
functions |14| | , which we will show can be used to compute 
the ALC. 

We introduce the correlators describing the fluctuation 
in Cooper pair occupation: 

Cj = yJ(n])-(n j f = y/(2-n j ) (nj). (6) 

These correlators can be directly evaluated for the vari- 
ational wavefunction @ giving Cj — 2\ujVj\ — Cj, so 
the local concurrence is precisely the fluctuation of the 
Cooper pair occupation. These fluctuations are localised 
within the range A of the Fermi level as given by 
The behaviour of the ALC is strongly influenced by the 
density of levels about the Fermi level, as depicted in 
Fig. 1. Hereafter we will only consider the case of uni- 
form distribution of the levels to simplify the analysis. 

In the thermodynamic limit @ becomes the exact 
ground state and the finite size corrections are of order 
1 /L ■ For systems with large but finite electron num- 
ber in the SC regime, the entanglement of the state © 
and the entanglement of the ground state will be the same 
up to corrections of order 1/L. However, the ground state 
cannot be accurately approximated by the product state 
(J3J in the FD regime. In this case the correlations spread 
out in ener gy sp ace over the entire width 2wd about the 
Fermi level [l4j . 

The canonical ensemble. Our next step is to establish 
that the correlators © are still equivalent to the local 
concurrences for a canonical system. Recall that since 
there are no unpaired electrons in the ground state, we 
can treat each ququadrit as an effective qubit. We ex- 
press the action of the canonical Fermi algebra on the 



FIG. 2: Ground state ALC for systems of L — 24 (dot-dash), 
40 (dash) and 68 (solid) levels. The bold and dotted lines 
show the analytic results for the thermodynamical limit and 
L — 2 respectively. The inset, showing the result for small A, 
highlights that the ALC is Independent and not universal. 



subspace of the Hilbert space with no unpaired electrons 
in terms of the Pauli matrices through the identification 
n = I + o z , c\_cl_ = cr + , c_c + = cr~. The uniqueness of 
the ground state and the u(l) invariance of the Hamil- 
tonian @ due to conservation of total electron number 

means that Q reduces to Cj = *Jl — (o- z ) 2 . Next we 
express the correlators (0 in terms of the Pauli matri- 
ces, with the result being Cj = Cj. A clarifying point 
is needed here. The local concurrence Cj is a measure 
of the entanglement between the effective qubit associ- 
ated with the jth level and the remainder of the system. 
Treating the jth level as a ququadrit, the reduced density 
matrix becomes 

Pi = ~(2 - (nj)) |00) (00| + ~ (nj) |11) (11| . (7) 

Taking the partial trace over either qubit yields a reduced 
density matrix which has the same non-zero eigenvalues 
as Q, so the local concurrence for either qubit within 
a ququadrit is the same as the local concurrence of the 
ququadrit viewed as an effective qubit. 

We see that the definition for the ALC in terms of the 
correlators © is the same for the canonical and grand 
canonical cases. Note that the derivation of the ALC in 
terms of © for canonical systems relied on invari- 
ance. In the thermodynamic limit the u(l) invariance of 
the ground state density matrix is broken, but the same 
expression for the ALC in terms of © is valid. 

To analyse the ground state ALC in the canonical case 
we employ results from the exact solution 0. An eigen- 
state of (j2J with M Cooper pairs is characterised by a 
set of complex parameters {v\ 1 ...,%} which provide a 
solution to the set of Bethe Ansatz Equations (BAE) 

2^1 ^2 

dX ^ Vl ~e k j£ v > ~ v * 
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io- 3 10* lnX io-' 10° 

FIG. 3: The ratio C /E versus ln(A). The results shown are 
for L — 24 (dot-dash), 40 (dash) and 68 (solid) levels, while 
the bold and dotted curves are the analytic results obtained 
for the thermodynamic limit and L — 2. The maximum for 
each case is a threshold coupling, below which other EMs 
become significant. For L = 2 and the thermodynamic limit 
the maximum occurs at A = 0, as there is only bi-partite 
entanglement in these cases. 

and the energy is given by E = 2^. =1 Vj + XdM. For 
the simple case of L = 2, M = 1, (JSJ) can be solved 
analytically yielding the ground state energy E = e± + 
€2 — \J <P\ 2 + (ei — £2) 2 - Using the Hellmann-Feynman 
theorem the correlators © can be computed which gives 
the ALC as C = A/VT + A 2 . 

For the case of general finite L the ALC can be com- 
puted through determinant representations of the Cj 
which are given as functions of the set . . . , vm}- The 
explicit formulae can be found in 0, |n| . fn the limits 
of weak and strong [2(j coupling for large but finite 
L we have the asymptotic results 

E ~ A/2(l + 0(l/L)), 

C ~ 1- 1/(6A 2 )(1 + 0(1/L)), A-^l, 

E ~ A 2 ln(2)/L(l + 0(l/lnL)), 

C ~ Xy/2/L ln(3 + v / 8)(l + 0(l/]ja.L)), A« 1. 

It is found that for A" 1 < 1, cf/E ~ 2/A, and 
C 2 /E ~ (21n 2 (3 + V8))/ln(2) « 8.97 



for A -C 1. Therefore in the FD and SC regimes the 
the quantity C /E displays scaling behaviour, i.e. the 
leading term is independent of L. In contrast to the FD 
and SC regimes, which are characterised by the scales 
l/ln(L) and l/L respectively, the scale l/ln 2 (L) occurs 
for the crossover regime 19]. 

In Fig. 2 we show the ALC versus A for L = 
2, 24, 40, 68 and the thermodynamic limit. Like the con- 
densation energy, it is clear that the ALC is a smooth 
monotonic function of A as it crosses from the FD to 
the SC regime. In Fig. 3 we plot the ratio C /E versus 
ln(A), for the finite cases L = 24,40,68. For sufficiently 
small A this ratio is approximately the constant value 8.97 
(up to a small correction 0(1/ ln(L))) in agreement with 
the asymptotic result, but is nonetheless monotonically 
increasing. At sufficiently large A the ratio is monotoni- 
cally decreasing and is well approximated by the analytic 
curve for the thermodynamic limit. For L — 2 and the 
thermodynamic limit the curves, also shown, are mono- 
tonic, while in each of the other three cases there is clearly 
a maximum at a finite value of ln(A). The coupling at 

which the maximum of C /E occurs is a threshold: below 
this coupling @ no longer approximates the ground state 
and we must expect that other (generally multi-partite) 
EMs become significant. 

For larger values of L, we appeal to a heuristic ar- 
gument based on the observation made in [l9l |: for the 
crossover regime, the condensation energy is roughly re- 
produced by simply summing up the contributions from 
the perturbative result in the FD regime and the BCS 
mean field theory in the SC regime |2l| . This also applies 
to the ALC. Therefore, the same picture we have drawn 
from exact Bethe ansatz solutions for small L works for 
very large L, thus filling in the gap between the tractable 
but relatively small Us and the thermodynamic limit. As 
L becomes very large, the threshold coupling tends to the 
value 2/ln(£) (the coupling at which mean-field theory 

breaks down 0). This gives C /E ~ ln 2 (X) at the 
threshold, showing the peak in Fig. 3 is not bounded as 
L increases. The competition of different scales in the 
crossover regime leads to the breakdown of the scaling 
behaviour observed in the FD and SC regimes. 
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